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O ! Abstract 
<D 

, Spectropolarimetry results for the starburst galaxy M82 are presented. The optical emission lines of the 

filaments in the energetic outflow ( "superwind" ) from the nuclear starburst region of M 82 are substantially 
polarized. The Ha polarization degrees and angles measured by our study are consistent with previous 
narrowband imaging polarimetry data. The polarized emission lines are redshifted with respect to the 
\ emission lines in the total light and systemic motion of the galaxy. The emission line intensity ratios 



[N n]/Ha and [S n]/Ha in the polarized light are similar to those of the nuclear star- forming region. In 
addition, the electron density N c derived from the [S n]A6731/A6717 line ratio of the polarized light is 
~ 600 — 1000 cm -3 at a distance of more than 1 kpc from the nucleus, whereas the N e derived from the 
Qh' total light are less than 300 cm' 3 . These facts strongly suggest that the emission from the nuclear starburst 

q , of M82 is scattered by dust grains entrained and transported outward by the superwind. A simple hollow 

biconical outflow model shows that the velocity of the outflowing dust grains, v<i, ranges from 100 to 200 
c/3 | km s _1 near the nucleus, decreases monotonically with the distance from the nucleus, and reaches ~ 10 

km s _1 at around 1 kpc. The motion of the dust is substantially slower than that of both ionized gas 
(«Hq ~ 600 km s _1 ) and molecular gas (i>co ~ 200 km s" 1 ) at the same distance from the nucleus of M82. 
This indicates that dust grains in the superwind are kinematically decoupled from both gas components 
' at large radii. Since the dust velocity v<± is much less than the escape velocity of M82 (w esc w 170 km s _1 

at 1.5 kpc from the nucleus), most of the dust entrained by the superwind cannot escape to intergalactic 
space, and may fall back into the galaxy disk without any additional acceleration mechanisms (such as 
radiation pressure) . 
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Introduction consumption of the interstellar medium by a superwind 

will cause star formation in a galaxy to cease abruptly, 
Active star formation in starburst galaxies plays one of although part of the expelled gas returns to the disk and 



the most important roles in galaxy evolution. In spite induces further star formation. These negative and pos- 

^ ■ of its short duration (~ 10 7-8 yrs), a starburst reforms itive feedbacks greatly affect the chemical evolution of a 

a substantial amount of the interstellar medium (ISM) galaxy. 

of a galaxy into stars. In addition, a starburst creates A significant amount of dust is associated with super- 
an enormous hot outflow from a galaxy via the collective winds. Submillimeter (sub-mm), far-infrared (IR), and 
effect of supcrnovae explosions and the stellar winds of mid-IR maps, optical color maps, and polarization stud- 
massive stars. This energetic outflow, known as a "super- ies of several nearby starburst galaxies have exhibited ex- 
wind," is a ubiquitous phenomenon in starburst galaxies tended dust emission along the galaxy minor axis (Alton 
(e.g., Heckman 2003; Veilleux et al. 2005). Theoretical et al. 1999; Leeuw & Robson 2009; Ichikawa et al. 1994; 
studies estimate that the terminal velocity of the hot gas Scarrott et al. 1991; Engelbracht et al. 2006; Kaneda et 
of a superwind reaches ~ 10 3 km s _1 (e.g Chevalier & al. 2010). Sub-mm observations of some starburst galax- 
Clegg 1985), far exceeding the typical escape velocity of ies suggest that the mass of the dust may reach ~ 10 6-7 
a galaxy. Metal-rich gas from the starburst region and M Q (Alton et al. 1999; Leeuw & Robson 2009). Although 
ambient disk gas are expelled from the galaxy disk and such a large-scale dust outflow has been inferred to play an 
pollute the galaxy halo and intergalactic space. Rapid important role in the evolution and metal enrichment of 



. the intergalactic medium and halo g well as the evo- 

Based on data collected with the Subaru Telescope operated by ... r , , .. lr r , r . , . , 

, , n , . c , lution of the host galaxy itscll, the late of this dust has 

the National Astronomical Observatory ot Japan .... 

thus far remained uncertain. This is because kinematic 
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information on dust outflow is quite difficult to obtain, 
since dust produces no sharp emission or absorption lines 
by which its radial velocity can be measured. 

One of the most promising techniques for probing the 
motion of dust in and around superwinds is optical spec- 
tropolarimetry. Outflowing dust grains entrained by a 
superwind scatter and polarize the continuum light and 
emission lines emanating from the nuclear starburst re- 
gion. In other words, the dust grains act as "moving mir- 
rors" for nuclear light. Hence, the velocity measured by 
the polarized emission lines of a superwind must reflect the 
motion of the dust with respect to the nucleus. Motivated 
by this idea, we carried out deep optical spectropolarime- 
try of the prototypical starburst galaxy M82 to reveal the 
dust kinematics of its superwind. 

M82 is distinguished by its very bright, kpc-scale su- 
perwind (Nakai et al. 1987; Seaquist & Clark 2001; 
Matsushita et al. 2005; Bierao et al. 2008; Bland & Tully 
1988; Shopbell & Bland-Hawthorn 1998; Ohyama et al. 
2002; Mutchler et al. 2007; Bregman et al. 1995; Tsuru et 
al. 2007; Ranalli et al. 2008; Strickland & Heckman 2007). 
The superwind of M82 is accompanied by large extrapla- 
nar dust filaments, as well as hot ionized gas (Ichikawa et 
al. 1994; Alton et al. 1999; Thuma et al. 2000). Alton et 
al. (1999) found a huge dust envelope extending along the 
minor axis of M82 using the 850- /im sub- mm observation. 
Recently, Leeuw & Robson (2009) detected a much fainter 
dust emission as far as 1.5 kpc from the galaxy disk. They 
found that the sub-mm morphology has a north-south 
asymmetry, which is consistent with the Ha and X-ray 
morphologies. Mid- and far-IR maps recently obtained 
using IR space telescopes reveal a complex structure in 
the kpc-scale filaments of polycyclic aromatic hydrocar- 
bon (PAH) dust extending along the minor axis of M82 
(Engclbracht et al. 2006; Kaneda et al. 2010; Roussel et 
al. 2010). The highly polarized nature of the optical con- 
tinuum and Ha emission of the outer region of M82 also 
indicates the presence of a vast quantity of dust in the 
superwind (Schmidt et al. 1976; Bingham et al. 1976; 
Visvanathan & Sandage 1972; Scarrott et al. 1991). The 
polarization degree of the Ha emission reaches 30% in 
some areas, and the polarization angle is almost perpen- 
dicular to the radial direction drawn from the nucleus of 
the galaxy (Scarrott et al. 1991). Therefore, M82 is an 
ideal object for studying dust kinematics in a starburst 
superwind via optical spectropolarimetry. 

In this paper, we present the results of deep optical 
spcctropolarimctric observations of the superwind of M82. 
This is the first attempt to ascertain the spatial structure 
of the dust motion in a starburst superwind via spectropo- 
larimetry. We adopted 3.89 Mpc as the distance to M82 
(Sakai & Madore 1999), which yields a linear scale of 18.9 
pc arcsec -1 for the galaxy. 

2. Observations 

Spcctropolarimctric observations of M82 were made 
with FOCAS (Kashikawa et al. 2002), attached to the 
Cassegrain focus of the Subaru Telescope (Kaifu et al. 




Fig. 1. The slit positions overlaid on a false color (bluc:V, 
green:_R, and red:Ha) image of M82. Regions discussed in the 
text are labeled R1-R8. 

2000), on December 22 and 23, 2003. The observations 
were carried out using the spectropolarimetric mode of 
FOCAS (Kawabata et al. 2003). We used a slit mask with 
eight 0".6 (width) x 20".6 (length) slitlets at 23". 7 inter- 
vals, and a VPH grism with 665 grooves mm -1 and a cen- 
ter wavelength of 6500 A. The resultant spectral resolving 
power was A/AA ps 1700, determined by the combination 
of the slit and the grating. The separation direction of 
the beam splitter was perpendicular to the direction of 
the slit length, and spectra of both ordinary and extraor- 
dinary rays were obtained simultaneously. 

The PA of the slit was set at 150° and 134°. We 
placed the slit center at the position of the 2.2 /im nu- 
cleus (Telesco et al. 1991). A unit data set consists of 
exposures taken at four different position angles (PAs; 0°, 
22.5°, 45°, 67.5°) of the half-wave plate. We took one data 
set with an exposure time of 600 s per exposure for PA 
150°, and four data sets with an exposure time of 720 s 
per exposure for PA 134°. Thus, the total exposure times 
were 2400 s for PA 150° and 11520 s for PA 134°. Figure 1 
shows the positions of the slits overlaid on an image taken 
with FOCAS (Ohyama et al. 2002). The regions in which 
the spectra were extracted and summed are also shown in 
Figure 1. 

We observed both an unpolarized star and a highly po- 
larized star to calibrate the polarization. We also observed 
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Fig. 2. Spectra of total light and polarized light around the Ho+[N II] wavelength region for R1-R4. The polarized spectra are 
represented by the solid line in the upper panels of each figure. The total light spectra are scaled by a factor shown in the upper-right 
corner of each figure and are represented by the dashed line. The polarization degree and polarization position angle are shown in 
the middle and bottom panels of each figure, respectively. 
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Fig. 3. Spectra around the Ha+[N n] wavelength region for R5-R8. The layout is the same as in Figure 2. 



a spectrophotometric standard star. The first night of the 
observing run, December 22, was clear and photometric. 
The second night was not photometric, with some small 
passing clouds. The seeing was approximately 0".6 the 
first night and 0".7~ 0".8 the second night. 

3. Data Reduction 

Standard CCD data reduction was applied to the raw 
CCD frames using IRAF. The bias level was subtracted 
from all frames. Since the dark current was negligible (< 1 
ADU), dark subtraction was not performed. Flat fielding 
was carried out using the dome flat data averaged over 
the four PAs of the half-wave plate. 

We extracted one-dimensional spectra of ordinary and 
extraordinary rays in eight regions (Rf-R8) of the two- 
dimensional spectral images. The locations of R1-R8 are 



shown in Figure 1. The distances from the nucleus and 
the areas of the regions are summarized in Table 1. We 
estimated the sky background components from the two- 
dimensional spectra in the outer region (140 - 145" from 
the nucleus), where no significant Ha flux was recorded. 
The polarization parameters were calculated from eight 
spectra (ordinary and extraordinary at the four wave- 
plate PAs) employing the method described in §6.1.2 of 
Tinbcrgcn (1996). Our observations of unpolarized stars 
indicated that instrumental polarization was negligible 
(< 0.1%). Moreover, our measurements of flat-field lamps 
through fully polarizing filters showed that the depolar- 
ization factor was also negligible 0.05). Therefore, we 
made no correction for instrumental polarization and de- 
polarization. The zero point of the position angle on the 
sky was determined from the observation of the highly 
polarized star. 
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Fig. 4. Spectra around the [S II] wavelength region for R1-R4. The layout is the same as in Figure 2. 



We fitted Gaussian profiles to the emission line spec- 
tra for both the total light and the polarized light. First 
we applied Gaussian fitting to the Ha+[N n] spectra, as- 
suming that Ha, [N il]A6548, and [N n]A6583 have the 
same radial velocities and FWHMs. We fixed the emis- 
sion line intensity ratio [N n] A6583/A6548 at 3 for this 
procedure. The Ha+[N n] emission lines of the total light 
in R3, R7, and R8 have a remarkable red-wing compo- 
nent, while those in R4, R5, and R6 exhibit a double- 
peak profile. We fitted double Gaussian profiles to the 
Ha+ [N n] lines of these regions to decompose their asym- 
metric profiles. The radial velocities and FWHMs of all 
components were assumed to be the same for the de- 
composition. We then applied Gaussian fitting to the 
[S ll]AA6717,6731 lines, assuming their radial velocities 
and FWHMs to be the same as those of the Ha+[N n] 
lines. Some asymmetric features were also found in the 



polarized emission line profiles of R3-R7, but the S/N ra- 
tios were not high enough to obtain a reliable decomposi- 
tion. Hence we applied a single Gaussian fitting to all the 
polarized emission lines, and then corrected for instrumen- 
tal broadening by applying the following simple equation: 



FWHM = y FWHMg t - FWHMf nst , where FWHM fit and 
FWHM inst are the FWHM obtained from the fitting and 
that of the instrumental profile, respectively We deter- 



mined that FWHM ir 



140 km s by fitting the sky 



emission lines near Ha. The physical parameters derived 
by the above procedure are listed in Tables 1 and 2. 

4. Properties of the polarized emission lines 

Figures 2 and 3 show the total fluxes, polarized 
fluxes, polarization degrees, and polarization angles of the 
Ha+[N ii] regions of the eight selected regions along the 
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Fig. 5. Spectra around the [S II] wavelength region for R5-R8. The layout is the same as in Figure 2. 



slits (R1-R8). Figures 4 and 5 show the same data for the 
[S n] region of the spectra. We summarize the physical 
parameters derived from the total light spectra in Tabic 1. 
Since double Gaussian fittings were applied to the emis- 
sion line profiles of R3-R8, the fitted parameters for each 
component are listed separately for these regions in Table 
1. The emission line fluxes of the polarized Ha, [Nn]A6581 
and [S n]AA6717,6731, the velocities and FWHMs of the 
polarized Ha, and the electron densities N c derived from 
the polarized [S n] line intensity ratios are given in Table 2. 
We assumed an electron temperature of 10 4 K in calculat- 
ing N c . The radial velocities are corrected to heliocentric 
values. 

The intensity- weighted polarization degrees of the emis- 
sion lines range from 5% to 15% (see Table 2). The po- 
larization degree begins to decrease on the blue side of 
the emission lines, reaches its minimum near the center of 



the lines, increases again on the red side of the lines, and 
attains a peak at the red wings of the lines. Except for 
Rl and R2, the polarization angles arc almost constant 
(« 40° - 60°) throughout the Ha+[N n] spectral region, 
and no difference exists between the polarization angles 
of the continuum and emission lines. In the spectra of 
Rl and R2, the polarization angles increase toward the 
red part of the emission lines (upper panels of Figure 2). 
The polarization angles are ks 35° in the continuum and in 
the blue part of the emission lines, while the angles reach 
~ 50° - 60° at the red wings of the lines. Except for the 
fine structures seen in Rl and R2, the polarization vectors 
are generally perpendicular to radial lines drawn from the 
M82 nucleus. 

Figure 6 shows the radial velocities of the total light 
and polarized light versus the distance from the nucleus. 
Clearly the polarized emission line spectra are redshifted 
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Table 1. Physical Parameters of the Ha-emitting Gas in the M82 Superwind 



ID 


PA* 


Dt 


Area 


hJ 


V 


FWHM 


/[N II] * 


f[S 1116717* 


f[S II] 6731* 


iVe,t 




(deg) 


(PC) 


(pcxpc) 




(km s- 1 ) 


(km s -1 ) 








(cm -3 ) 


Rl 


150 


140 


11x28 


97 


215 ±10 


189 ±22 


46 


15 


13 


370 ±10 


R2 


134 


245 


11x32 


190 


234 ±10 


161 ±26 


85 


26 


25 


450 ±10 


R3 


150 


630 


11x38 


7.5 


192± 11 


150 ±28 


3.3 


1.4 


1.1 


120 ±10 










2.3 


385 ±12 


224 ±19 


1.4 


0.4 


0.4 


370 ±10 


R4 


150 


735 


11x28 


6.7 


195 ±11 


187 ±23 


2.9 


1.1 


0.85 


110±10 










2.6 


-32 ±14 


75 ±63 


0.83 


0.41 


0.31 


120 ±10 


R5 


150 


805 


11x32 


8.7 


205 ±11 


121 ±35 


4.4 


2.1 


1.6 


120 ±10 










2.4 


-43 ±13 


95 ±46 


0.67 


0.31 


0.26 


260 ±10 


R6 


150 


910 


11x43 


7.7 


215± 11 


116 ±37 


3.6 


1.7 


1.2 


70 ±10 










1.7 


-40 ±13 


107 ±41 


0.57 


0.24 


0.21 


320 ±10 


R7 


134 


840 


11x29 


3.0 


83 ±10 


178 ±23 


1.2 


0.75 


0.50 


< 50 










2.5 


167± 11 


325 ±12 


1.7 


0.41 


0.37 


400 ±10 


R8 


134 


980 


11x77 


3.0 


75 ±14 


130 ±33 


1.3 


0.98 


0.61 


< 50 










7.0 


155 ±11 


295 ±14 


4.3 


1.4 


1.2 


260 ±10 



* The position angle of the slit. 

t The projected distance from the nucleus. 

t The emission line flux in units of 10 — 15 erg s -1 cm -2 . 



Table 2. Physical Parameters of the Gas Measured in Polarized Emission Lines of the M82 Superwind 



ID 


P.Deg.* 

(%) 


P.PAT 

(deg) 




Up 

(km s- 1 ) 


FWHM 

(km s -1 ) 


/[N II]6583,p* 


f[S II]6717,p* 


f[S II]6731,p* 


^c,p 

(cm- 3 ) 


Rl 


4.5 


35 


44 


252 ±11 


219 ±19 


24 


4.2 


4.8 


890 ± 10 


R2 


3.4 


37 


65 


247 ±11 


203 ±21 


33 


6.7 


7.6 


910 ±10 


R3 


12.7 


65 


12 


317 ±14 


279 ±15 


6.4 


0.70 


0.67 


520 ±50 


R4 


10.5 


49 


9.8 


240 ± 15 


268 ±16 


4.6 


0.51 


0.40 


130 ±50 


R5 


11.4 


51 


12 


226 ±13 


221 ±19 


6.9 


1.2 


1.2 


740 ± 40 


R6 


9.8 


54 


8.6 


227 ±15 


221 ±19 


4.1 


0.44 


0.47 


750 ±130 


R7 


10.4 


38 


5.6 


213 ±14 


312 ±14 


3.0 


0.50 


0.44 


310 ±40 


R8 


15.4 


36 


15 


234 ±13 


240 ±18 


8.1 


1.3 


1.5 


940 ± 20 



* The intensity-weighted mean polarization degree of the Ho? emission. 

t The intensity-weighted mean position angle of the polarization vector of the Ha emission, 
t The polarized emission line flux in units of 10 -lfi erg s -1 cnr 2 . 



with respect to the total spectra (see also Figures 2-5). 
Also, the polarized Ha lines have radial velocities higher 
than the systemic velocity of M82 (v sys = 203 km s -1 ; 
Shopbell & Bland-Hawthorn 1998; Gotz et al. 1990). The 
velocity excess, v p — v sys , ranges from ~ 10 km s -1 to ~ 100 
km s -1 . In the regions close to the nucleus (Rl and R2) 
v p ~ v sya ~ 50 km s , increasing to ~ 100 km s -1 in R3. 
Note that R3 is in the northwest part of the superwind of 
M82. In the southeast regions (R4-R8), whose distances 
from the nucleus are 0.7 - 1 kpc, the polarized Ha lines 
have almost constant velocities, and v p — v sys ~ 10 km s . 
In contrast, the Ha lines of the total light are blueshiftcd 
relative to v sys in R4-R8, which is consistent with previous 
studies (Shopbell k Bland-Hawthorn 1998; Greve 2004). 

The emission line intensity ratios [N n]/Ha and 
[S n](A6717±A6731)/Ha for the polarized light are similar 
to those of the nuclear region. Figure 7 shows a plot of the 
ratios for the total light and the polarized light, together 
with the data taken around the nucleus (Westmoquette et 
al. 2009). The data for the polarized light are distributed 
in and around the range of the nuclear data. The data for 



the total light, however, are distributed well outside the 
nuclear data (Figure 7). 

The emission line intensity ratios [S n]A6731/ 
[S n]A6717 of the polarized spectra are higher than those 
of the total light in all regions except R4 (Figures 4 and 
5). Figure 8 shows the N e derived from the [S n] line ra- 
tios. The electron densities iV e ,t obtained from the [S n] 
line ratios of the total light are ps 400 cm -3 at ~ 0.2 kpc 
from the nucleus and decrease to < 300 cm -3 at >0.5 
kpc (Figure 8). In contrast, the electron densities N e>p 
measured in the polarized spectra arc almost constant (~ 
500 - 1000 cm -3 ), except in the two low-density regions 
(R4 and R7). This indicates that scattered light origi- 
nates from a high-density region. The electron density of 
the nucleus of M82 is « 1000 cm -3 (Westmoquette et al. 
2009). The iV 0i p we derived are slightly smaller than, but 
almost consistent with, the nuclear N e . 

The foregoing polarized light characteristics (polariza- 
tion angles, emission line ratios, and N e , p measured in the 
polarized light) are consistent with the idea that the po- 
larized light is nuclear light scattered by the dust in the 
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Fig. 6. Position-velocity diagram of the Ha emission of the 
superwind of M82. The filled red and blue circles indicate the 
heliocentric radial velocities of the polarized Ha emission and 
the intensity-averaged velocities of the Ha emission for the 
total light, respectively. The profiles of the total light emis- 
sion lines in the outer regions (R3-R8) are fitted with double 
Gaussian functions, and the velocities of the split components 
are then plotted as open blue circles. The dotted line repre- 
sents the systemic velocity of M82. 
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Fig. 7. The emission line intensity ratios of the polar- 
ized light and the total light [S il] (6717+6731) /Ha vs. 
[N Il]6584/Ha. The filled red and blue circles indicate the 
data for the polarized light and the total light, respectively. 
The open blue circles represent the data for the split compo- 
nents of the total light. The open squares show the data in an 
area of Ri 50 X 100 pc around the nucleus of M82, measured by 
integral field spectroscopy (Westmoquette et al. 2009). The 
upper and lower squares represent the data for the narrow 
component and the broad component, respectively ("CI" and 
"C2" in Westmoquette et al. (2009)). The error bars associ- 
ated with the squares show the distribution of the data over 
the measured region. 



Fig. 8. Electron density derived from the [S II] emission lines 
of M82. The filled red and blue circles indicate the data for 
the polarized light and total light, respectively. The open 
blue circles represent the data for the split components of the 
total light. The dashed line marks the lower limit of electron 
density that can be derived from the intensity ratios of the 
[S II] emission lines. 

superwind of M82. In addition, we discovered that the 
radial velocities of the polarized emission lines are sys- 
tematically greater than the velocity v sys of the galaxy. 
The dust grains that scatter the nuclear light can be in- 
ferred to move outward from the nucleus. Accordingly, 
our research strongly suggests that dust grains entrained 
by the superwind arc transported outward from the disk 
by the wind and scatter the nuclear light. 

5. Outflow velocity of the dust wind 

We attempt to derive the true outflow velocity of the 
dust grains in the superwind of M82 using the polarized 
spectra. Assuming that the radial velocities of the polar- 
ized emission lines reflect the motion of the dust grains 
with respect to the nucleus, we construct a hollow biconi- 
cal outflow model for the dust outflow of M82 (Figure 9). 
Since forward scattering is much more efficient for scat- 
tering by ordinary interstellar dust grains (van dc Hulst 
1957), we ignored the backscattering component from the 
far sides of the outflow cones. 

The idea of the model is as follows. Dust grains are 
flowing outward from the galaxy disk along the walls of 
the hollow cones, with outflow velocity va- The radius of 
the base of the outflow cone is b kpc. The dust grains act 
as moving mirrors, scattering the light from the nucleus. 
The opening angle of the outflow is 2x9, and the angle of 
inclination of the axis of the outflow cone is i (Figure 9). 
Although i is generally different from the angle of inclina- 
tion of the galaxy disk (i.e., the axis of the cone may be 
tilted with respect to the rotational axis of the galaxy) , the 
difference between the two values should be small, as sug- 
gested by the ionized gas wind morphology (Shopbell & 
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Bland-Hawthorn 1998; Ohyama et al. 2002; Greve 2004). 
Hence we consider i to be effectively the same as the incli- 
nation of the galaxy. In the following calculation, the sign 
of i is negative on the northwestern side of the galaxy and 
positive on the southeastern side. When the dust grains 
are observed at a projected height h above the nucleus, 
the distance z between the observed dust grains and the 
galaxy disk ("true height") is expressed as follows: 

h + b tan(i) cos(i) 
2 ~~ 1 - tan(0) tan(i) ' ^ ' 

The distance r<\ between the position of the dust grains 
projected onto the galaxy disk and the nucleus (i.e., 
the radial distance within the galaxy disk) is then = 
z tan(f?) + b. Using these values, we can express the angle 
ip between the outflow velocity Vd and the line connecting 
the observed flow region and the nucleus as 



(2) 



axis of the dust 
outflow cone 
/ 



%j) = 7i72- 8 -tan -1 — 
Jd 



The dust grains are receding from the nucleus with a ve- 
locity I'd cos(V0 (red arrows in Figure 9), and the scattered 
emission lines are redshifted with this velocity. In con- 
trast, the velocity component of the grains along the line 
of sight is Vd sm(8 + i) (blue arrows in Figure 9), and hence 
the scattered emission lines, are blucshiftcd with this ve- 
locity, relative to the systemic motion of the galaxy. As 
a result, the observed radial velocity v p of the scattered 
emission lines will be v p — Vd cos(ip) — {vd sm(8 + i) — v sys }. 
In this simple model, we can thus derive the dust outflow 
velocity Vd as follows: 



cos(V ? ) — sin(0 + i) ' 



(3) 



where v p and v sys (= 203 km s" 1 ; Gotz et al. 1990) are 
the radial velocity of the polarized emission lines and the 
systemic velocity of M82, respectively. 

To calculate Vd, we make the following assumptions 
about the parameters, b, 8, and i. The dust can be rea- 
sonably expected to be entrained from the galaxy disk 
outside the main hot outflow. The radius of the base of 
the dust outflow is then assumed to be the same as the 
radius of the central starburst, so that b ps 0.3 kpc (Greve 
2004). The opening angle 28 and the angle of inclination 
i of the dust outflow are hardly determined only by our 
observations. We examine two cases for 9: a narrow-angle 
case and a wide-angle case. In the narrow-angle case, we 
assume that the dust grains move in the same direction as 
the ionized gas. In this case, we adopt the value 20 = 25°, 
which was derived by Greve (2004). If the dust flow is 
associated with molecular gas, 8 may be much wider than 
the ionized gas flow (the wide-angle case) . When the dust 
grains are associated with the molecular outflow identified 
by Walter et al. (2002), 28 « 50°. Ohyama et al. (2002) 
suggested a wider opening angle for the dust outflow in 
or near the galaxy disk (28 > 90°) using an optical extinc- 
tion map (see also Ichikawa et al. 1994). However, the 
dust scattering the nuclear light is very far away from the 
disk according to the imaging polarimctry maps (Scarrott 
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Fig. 9. A schematic diagram of the M82 dust wind. The 
dust grains in the galaxy disk are assumed to be entrained by 
the superwind expelled from a circumnuclcar circular region 
whose radius in the disk is b. The entrained dust flows along 
the walls of a cone whose opening angle is 2 X 6. The angle 
of inclination of the axis of the dust flow cone is i. The dust 
grains in the wind reflect the nuclear emission, acting as mir- 
rors moving at velocity v^cosip with respect to the nucleus, 
where is the outflow velocity of the dust. 

et al. 1991). Hence the wider-angle dust flows near the 
disk plane would not contribute to the polarized light we 
detected here. Thus, we adopt the value 28 = 50° in the 
wide-angle case. Finally, wc assume that i = 15°, which 
was derived by Greve (2004), for both cases. 

Wc find that the estimated dust outflow velocities based 
on the above simple model are much less than the escape 
velocity of M82. In Figure 10, the dust outflow velocities 
«d of all regions arc plotted against the height z measured 
from the galaxy disk for the above two cases. The dust 
outflow velocity Vd decreases almost monotonically with 
z and reaches a minimum value of ~ 20 - 30 km s _1 
at approximately 1-1.2 kpc in both cases. Martin (1998) 
estimated the escape velocity of M82 using the H I velocity 
field and found that v esc « 170 km s _1 at 1.5 kpc from 
the nucleus (Martin 1998). The present estimate for Vd is 
much less than this value. Some indication of an upturn in 
Vd for z > 1.2 kpc is observed (see Figure 10), but whether 
this trend reflects a global acceleration of the dust grains 
at high altitudes above the disk is not clear, owing to the 
lack of data beyond z ~ 1.5 kpc. 

We fit a constant deceleration model to the dust ve- 
locity field in the M82 superwind and estimate the dust 
recycling timescalc. Assuming that the dust grains are 
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Fig. 10. Estimated dust outflow velocity ty of M82. The 
filled squares and open squares represent the data derived 
for the narrow-angle case and the wide-angle case, respec- 
tively. The horizontal axis denotes the height measured from 
the galaxy disk (z in Figure 9). The solid and dotted lines 
represent best-fit curves of constant deceleration models (see 
text) for the narrow-angle case and the wide-angle case, re- 
spectively. 

entrained by the superwind at an initial velocity v^o 
and are decelerated at a constant rate d r , we fit linear 
functions Vd{z) — t>d,o + d r z for the two cases and find 
(vd,o km s -1 , d r km s^ 1 kpc -1 ) = (186, —139) for the 
narrow-angle case and (268, —176) for the wide-angle 
case. We can then calculate the distance z(t) (in kpc) 
attained by the dust grains at time t (in Myr) from 
z(t) = v ( i,o/d I -{l — exp(— 1.1 x 10~ 3 -d r i)}, since the grains 
are entrained from the galaxy disk. In this simple model, 
the dust grains reach a height of 1 kpc above the disk in 
~ 10 7 yr in both cases. If the deceleration rate d r is con- 
stant, «d would equal zero at ~ 1.5 kpc above the disk. 
The free-fall timescale from 1.5 kpc above the disk of M82 
to the disk plane is calculated as 2 - 3xl0 7 yr, under the 
assumption that a mass of ~ 10 10 M© is confined within 
2 kpc from the M82 nucleus (Sofue, Y. 1992). Hence, if 
the deceleration of the dust flow is monotonic, the dust 
expelled by the superwind may return to the disk in sev- 
eral times 10 7 yr. Note that this timescale is comparable 
to that of starburst activity. 

6. Radial extension of the dust outflow of M82 

Our simple dust outflow model suggests that dust grains 
cannot go beyond 1.5 kpc above the disk of M82. However, 
a body of evidence indicates that dust actually distributes 
over a few kpc around M82. The imaging polarimetry ob- 
servations of Scarrott et al. (1991) showed that the Ha 
emission is highly polarized even at the edge of the ion- 
ized gas flow of M82 (Scarrott et al. 1991). At the tip 
of the Ha emission region, ~ 3 kpc from the disk, the 



polarization degree reaches 30%, which means that a sig- 
nificant amount of dust exists at that altitude. Sub-mm 
observations have detected dust emission extending more 
than a few kpc from the disk (Leeuw & Robson 2009). 
An 8-/xm image taken with the Spitzer infrared satellite 
revealed a complex filamentary dust structure extending 
3 kpc from the disk (Engelbracht et al. 2006). Recently, 
very extended extraplanar cold dust emission in M82 was 
discovered via mid- and far-IR observations (Kaneda et 
al. 2010; Roussel et al. 2010). The UV emission extend- 
ing to 6 kpc along the minor axis of M82 has also been 
attributed to dust scattering (Hoopes 2005). 

How dust grains are transported to such a high altitude 
remains a mystery. If the dust grains originate from the 
galaxy disk, additional mechanisms for their outward ac- 
celeration (such as radiation pressure from the starburst 
or further mechanical acceleration by the wind) are needed 
to explain such an extended dust distribution. Fcrrara et 
al. (1991) calculated the effect of radiation pressure on the 
motion of interstellar dust grains. They determined that 
dust grains are expelled to intergalactic space at several 
hundred km s _1 by the radiation pressure of stellar light in 
our Galaxy and the star-forming galaxy NGC 3198. Their 
calculation predicted that the velocity of the dust would 
rapidly increase to ~ 500 km s _1 within 10 Myr from the 
onset of the flow and then decrease to ~100 - 200 km s _1 
in 100 Myr. The dust grains arc transported outward to 
an order of 100 kpc in 100 Myr. 

The results of Ferrara et al.'s calculation (Ferrara et 
al. 1991) clearly contradict our observational results for 
M82. Our observations suggest that Vd decreases with 
time after the onset of the flow (see Section 5). Note that 
radiation pressure is highly dependent on the dust grain 
size a. Radiation pressure is proportional to a 2 , whereas 
gravitational force is proportional to a 3 . Although the 
grain size distribution and spatial distribution of the ex- 
traplanar dust of M82 are not known, dust components 
traced by mid-IR observations (Engelbracht et al. 2006; 
Bierao et al. 2008; Kaneda et al. 2010) should consist pri- 
marily of small PAH dust particles. These small dust 
particles may be transported very far away from the disk 
by radiation pressure and form a highly protracted struc- 
ture. Dust particles involved in optical scattering, how- 
ever, are much larger and heavier than PAH dust; sub- 
micron-sized grains are mostly responsible for optical scat- 
tering (2na/X^ 1; e.g., van de Hulst 1957). Thus, the dis- 
crepancy between mid-IR dust extension and optical scat- 
tering dust kinematics is possibly reconciled by the size 
difference between the grains responsible for each emis- 
sion. 

Note, however, that the Ha emission is scattered even at 
the tip of the bright Ha filaments (Scarrott et al. 1991), 
indicating that optical scattering dust is still present at 
this distance. At least two possible explanations exist for 
this finding: 1) The origin of the outer dust is the same as 
that of the inner dust, and the grains are accelerated by 
some mcchanism(s) at > 1 kpc. 2) The origin of the outer 
dust is different from that of the inner dust. Some indi- 
cation of dust acceleration is observed at around 1.2 kpc 
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from the disk, since a velocity upturn appears in the veloc- 
ity field we obtained (see the previous section). This may 
indicate that the dust grains are accelerated outward be- 
yond 1 kpc. However, the velocity upturn is not especially 
significant, since the lack of data points makes it impos- 
sible to be certain whether the dust is really accelerated 
beyond 1 kpc. However, the distant dust may be associ- 
ated with the H I stream traversing the disk of M82 (Yun 
et al. 1993). This stream is part of a giant H I cloud com- 
plex extending across the M81/M82 group, and possibly 
formed by a galaxy-galaxy interaction between M82 and 
M81 (Applcton, Davies & Stephenson 1981; Chynowcth 
ct al. 2008). Recently, Roussel et al. (2010) detected a 
protracted cold dust emission around M82 via far-IR (250 
/jm, 350 /xm and 500 fim) observation with the Herschel 
satellite. The cold dust extends to 9 kpc from the center. 
By comparing their far-IR maps with other wavelength 
data, they found that the outer part of the cold dust is 
associated with the H I stream, while the inner part is 
dominated by the warm (~ 30 — 50 K) component asso- 
ciated with the superwind (Roussel et al. 2010). In this 
case, optical scattering dust far away from the disk is pro- 
vided by this stream, and the dust velocity field in the 
outer part of the superwind may reflect H I stream kine- 
matics. To distinguish between the two possibilities listed 
above, one must trace the motion of the dust toward a 
higher altitude above the disk. 

7. Comparison with other outflow materials 

We compare the kinematics of the dust and other out- 
flowing materials, including ionized gas, molecular gas, 
and neutral atomic gas, in the superwind of M82. 

The kincmatical dissociation between the dust and the 
warm ionized gas (optical emission line gas) of M82 is very 
clear. The outflow velocity of the ionized gas increases 
almost linearly with the distance from the disk at a rate 
of 300 - 400 km s" 1 kpc" 1 (Greve 2004). In contrast, the 
dust outflow decelerates at a rate of —150 km s -1 kpc -1 
(see Section 5). The final outflow velocity of the ionized 
gas in M82 reaches ^500 - 600 km s _1 at ~ 1 kpc from 
the disk (Shopbell & Bland-Hawthorn 1998; Greve 2004), 
which is much higher than the dust velocity we derived. 

The velocity field of the molecular gas in the M82 super- 
wind is also inconsistent with the kinematics of the dust. 
Several CO observations have detected a molecular gas 
outflow with velocity ~ 200 km s _1 in M82 (Seaquist & 
Clark 2001; Walter ct al. 2002). This value is almost the 
same as the dust velocity near the nucleus (z < 0.4 kpc) 
for both the narrow- and wide-angle cases of our dust out- 
flow model. Hence dust and CO gas may be associated 
with each other near the nucleus. The spatial coincidence 
between CO emission and the 450-^m emission also sup- 
ports a spatial association of dust and molecular gas in 
the central region (Hughes ct al. 1994). However, the CO 
outflow velocities (~ 200 km s _1 ) were measured at ~ 1 
kpc from the nucleus, where the dust velocity decreases to 
less than 100 km s _1 (Figure 10). Thus, we can conclude 
that molecular gas is also kinematically decoupled from 



the dust at high altitudes in the superwind of M82. 

Finally, we compare dust motion with neutral atomic 
gas motion. Hcckman ct al. (2000) discovered that Na I 
absorption lines are systematically blueshifted with re- 
spect to the systemic velocities of starburst galaxies (see 
also Martin et al. 2002, 2005; Rupke et al. 2002, 2005; 
Schwartz & Martin 2004). They interpreted this phe- 
nomenon as a neutral gas outflow from the starburst re- 
gion. Using a correlation between the optical absorption 
coefficient Ay and the equivalent width of the Na I ab- 
sorption line, they inferred that dust is associated with 
neutral atomic gas, which is represented by Na I absorp- 
tion. The outflow velocity derived from the blueshift of 
the Na I lines of M82 is ps 80 km s _1 . This value is signifi- 
cantly smaller than near the nucleus, which may imply 
that neutral atomic gas is also kinematically decoupled 
from the dust near the nucleus. However, we should note 
that Hcckman et al. (2000) measured the line-of-sight ve- 
locity of Na I only in the nuclear region. If the outflow 
were spherically symmetric, the Na I velocity they derived 
would equal the actual outflow velocity. However, various 
observations indicate that the superwind of M82 has a 
bipolar nature. Hence the Na I velocity represents not the 
extraplanar outflow velocity, but the expansion velocity 
of the circumnuclear gas in the galactic plane. In other 
words, currently available data do not allow us to directly 
compare the kinematics of dust and neutral atomic gas. 
Deep spectroscopy of Na I absorption lines at large radii 
in M82 would help us to study the detailed outflow kine- 
matics of neutral atomic gas and to compare it with dust 
kinematics. 

8. Conclusion 

We conducted optical spectropolarimetry observations 
of the starburst superwind of the prototypical starburst 
galaxy M82 to reveal the motion of the dust entrained by 
the superwind. The Ha polarization degrees (~5 - 15%) 
and angles measured by our study are consistent with pre- 
vious narrowband imaging polarimetry maps. The polar- 
ized emission lines are redshiftcd relative to the systemic 
motion of the galaxy. The emission line intensity ratios 
[N n]/Ha and [S n]/Ha in the polarized light are similar 
to those of the nuclear star-forming region. The electron 
densities N e derived from the polarized [S n] line ratio 
are much higher than those derived from the total light. 
These facts strongly suggest that the emission from a nu- 
clear starburst is scattered by dust grains entrained and 
then transported outward by the starburst superwind. We 
derived the outflow velocity of the dust grains, i>d, using 
a simple hollow biconical outflow model. The outflow ve- 
locity «d is on the order of a few hundred km s _1 near 
the nucleus and decreases monotonically with the distance 
from the nucleus. The dust motion revealed by this study 
is substantially slower than the motion of the other com- 
ponents of the superwind (ionized gas and molecular gas). 
The outflow velocity of the dust is also much less than the 
escape velocity of M82. In the absence of any additional 
effective acceleration mechanisms (such as radiation pres- 
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sure), the dust expelled by the superwind would fall back 
into the galaxy disk within several times 10 7 yr. 
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